Structural glass on a lattice in the limit of infinite dimensions. 
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We construct a mean field theory for the lattice model of a structural glass and solve it using the 
replica method and one step replica symmetry breaking ansatz; this theory becomes exact in the 
limit of infinite dimensions. Analyzing stability of this solution we conclude that the metastable 
states remain uncorrelated in a finite temperature range below the transition, but become correlated 
at sufficiently low temperature. We find dynamic and thermodynamic transition temperatures 
as functions of the density and construct a full thermodynamic description of a typical physical 
process in which the system gets trapped in one metastable state when cooled below vitrification 
temperature. We find that for such physical process the entropy and pressure at the glass transition 
are continuous across the transition while their temperature derivatives have jumps. 
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Thermodynamics of structural glasses and of vitrifica- 
tion transition is a long standing problem in statistical 
physics. The essential feature of glass formation is a divi- 
sion of the phase space into an exponentially large num- 
ber of similar compartments. The system gets trapped 
in one of these compartments and stays there forever in 
a complete defiance of ergodic hypothesis; this defines 
dynamical "states" of the system. Each of these states 
is characterized by a local order parameter (e.g. atoms 
positions, magnetization) that varies in space and distin- 
guishes one state from another. This makes even a mean 
field theory difficult to construct because the thermody- 
namic ground state of the system (or any other particular 
state) has a unique configuration of the order parameter 
field, thus making it difficult to describe system in terms 
of average (and site independent) quantities. Few ways to 
resolve this difficulty were proposed recently [[0-0; they 
all share a common idea that states of the system with 
the same energy per site (but perhaps different total en- 
ergies) are essentially equivalent and instead of studying 
one particular state one can average over all states with 
the same energy density. 

There are a few ways to implement averaging over 
states, the most convenient ones seem to be the clone 
method |5|,|6| and the introduction of a small random field 
conjugated to the order parameter (magnetization in case 
of spin glasses or density in structural glasses) 0. The 
physical idea of the latter is that an infinitesimally small 
random field applied to a system with many metastable 
states rearranges the energies of low-lying states mak- 
ing the problem similar to one with quenched disorder. 
In a spin system, for example, we add to the Hamil- 
tonian a magnetic field part: TL — > TL + J^ hiSi with 
small random hi. The resulting change in the energy 
of a typical metastable state is of the order of \f~Nh; 
because this energy interval contains a large number of 
metastable states, we expect that a small non-zero field 
would result in a large rearrangement of their energies 
but would not change the properties of individual states. 



Averaging over the random field configurations may be 
then performed in a usual way introducing n replicas of 
the system and taking the limit n — » 0. The alterna- 
tive idea of the clone method is to consider a system 
of m clones constrained to be in the same metastable 
states by an infinitesimally weak coupling. Generally, 
the free energy density of the clone system may be 
written as F clone = -TJV" 1 In/ dFe- N( - m()F - s °°»t ( - F '» 
where S con f{F) is configurational entropy (S con f — 
j; ln(Af states )). If S con f (F) is concave and dS conf (F)/dF 
is finite at the lower bound F — Fo (corresponding to 
the ground state) the appropriate choice of m (m oc T 
at low T) leads to a partition function dominated by a 
small vicinity of any given F > F$ that still contains 
thermodynamic number of states providing the averag- 
ing mechanism in this approach. In this paper we shall 
apply the mean field formalism developed in our paper 
M that combines random field approach and the locator 
expansion developed for the glass physics in |1 to the 
simplest model of a realistic structural glass. The mean 
field theory for the structural glasses can be justified for- 
mally only in the limit of infinite dimensions; thus, in our 
derivation of the mean field theory below we shall con- 
sider the space of arbitrary dimension, d, and keep only 
the leading terms in 1/d. We hope, however, that it pro- 
vides reasonable results even for d = 3 because the actual 
parameter is the coordination number which is already 
large in three dimensions. 

Another important physical assumption that simplifies 
the calculations significantly is that different states are 
uncorrelated. Formally, this allows one to look for one 
step replica symmetry breaking (RSB) solution of the 
mean field equations. We check the stability of this solu- 
tion below and find that it is stable in some temperature 
range below the glass transition temperature but that it 
becomes unstable at low temperatures. The results of one 
step RSB can be translated much easier into the physi- 
cal properties than results of a full RSB because in this 
case one can establish a formal equivalence between this 



method and a clone approach with the number of clones 
being equal to the size of the blocks in 1RSB ansatz m, 
and replica free energy F rep i being equal to the free en- 
ergy of the clone system per one clone F rep i = F c i one /m. 
For a given m the main contribution to F c i one comes from 
states with m = TdS con f(F)/dF , thus 
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Sconf = f3m 2 -^ :: (F clone /m), F = -^F cione . (1) 
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Metastable states first appear at the dynamical tran- 
sition temperature T g , at which the system becomes 
trapped in one of the metastable states with largest 
possible free energy F because these states dominate 
the configuration space since S con f{F) is a monotoni- 
cally increasing function. Under a further decrease of 
temperature the system remains trapped in the same 
metastable state, thus the thermodynamics of a physi- 
cal cooling process is determined by the free energy of a 
single metastable state. To describe this physical pro- 
cess by a replica method we need to find the depen- 
dence m(T) consistent with the requirement that the 
physical system is trapped in a single metastable state. 
We assume that all states that were similar at one tem- 
perature remain similar as the temperature is decreased 
and these states do not bifurcate or disappear (consis- 
tent with one step RSB); this implies that configuration 
entropy S con f is constant along the trajectory m(T) cor- 
responding to a physical cooling process. We check that 
this assumption is consistent with (fil) that gives the free 
energy of a metastable state and the corresponding con- 
figuration entropy in terms of the replica free energy. 
First, we note that the usual thermodynamic identity 
E = d[f3F rep i(f3,m)]/d(3 that relates internal energy to 
the free energy remains valid in replica approach. Sec- 
ond, if the system is trapped in a single metastable state, 
the thermodynamic identity E(/3,m) — d[f]F(f3,m)]/d(3 
must also hold. Remarkably, these equations are con- 
sistent with each other and with (|lj) if and only if 
SconfiP, m ) — const. This gives us the implicit equation 
for m(T) and allows us to deduce the thermodynamics of 
the system trapped in a single metastable state provided 
one step RSB is a stable solution. 

We now turn to the details of the model and its solu- 
tion. We consider the simplest glass model which con- 
tains only short range repulsion between atoms: 
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where i,j are the sites of a d-dimensional lattice and 
Pi = 1,0 represents the on-site density. The glass is 
formed for such coupling matrices Jij that do not allow 
low energy crystal states (density waves). It is conve- 
nient to represent Jij in the form J = uf{t/\/2d), where 
t is the nearest neighbor hopping operator and we intro- 
duced the scaling \/2d that is needed to get a sensible 



limit at d — > oo |Hj. The simplest choice of the func- 
tion /(e) ~ e leads to the crystal ground state with a 
chess-board ordering for density p rj 1/2. We shall take 
/(e) = e 2 and show that this frustrated interaction leads 
to the formation of the glass state at low temperature. 

The Hamiltonian (0) may be represented as H = 
Hq + H' , separating out large but physically irrelevant 
constant term Hq = udp 2 . We focus on the nontrivial 
term H' = J2i j P~iJi,jPj where p = p — p. Introducing 
the replicas we get 
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where a is the replica index taking values form to n, and 
p, a are the the Lagrange multipliers (equal to the chem- 
ical potential divided by the temperature at the saddle 
point) that impose constraint on the number of parti- 
cles in each replica. In the limit of infinite dimensions 
the original model (||) can be replaced with an effective 
single site model 



P H 'eff = o U ( B ) ~ ^PaBajpp + p a p a 
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where B is the matrix in replica space and p = p — p 
has only replica index. The mean field self-consistency 
condition is that B gives the saddle point of the free en- 
ergy: (p a pp) — dU{B)/dB a p. The potential U(B) can 
be determined from the condition that all single site cor- 
relation functions of the model (0) coincide with the cor- 
relation functions of the original model d|). Instead of 
comparing these correlation functions directly we decou- 
ple the interaction term introducing the auxiliary fields, 
v a , conjugated to p 
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and compare the correlation functions of these fields. 
Summation over p in both models results in the on-site 
interaction of the fields v. Inspecting the terms of the per- 
turbation theory in this on-site interaction for the corre- 
lator G"? = {vfVj) one verifies that in the leading order 

in 1/d it is given by G — [T J -1 — S] _1 with the self energy 
£ which is diagonal in the site index: £ = — {A(3u)^ 1 8ij. 
This approach is similar to a locator expansion 0] but in 
our case the locator A is non-trivial in the replica space. 
The single site correlation function Q a p = G°f (that we 
need to establish the correspondence between the mod- 
els) is 



g = u(l deu{e)[l/f(s) + A- 



(7) 



where v(e) — exp[— e 2 /2]/\/27r is the density of states in 
the limit d — > oo . 

Now we turn to the model (g). Here the self energy 
is diagonal in the site index by construction, further, the 
interaction part of this model is the same as for model 
(||); assuming that their single site correlation functions 
coincide, we get that their single site self-energies are 
equal, so that for this model 
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g = [B- 1 + (PuA)- 1 ]- 



(8) 



which gives us equation relating B and A. Further, the 
correlation function of the auxiliary fields can be ex- 
pressed through V = {papp) via Q = B — BVB, combin- 
ing this equation with the saddle point condition and 
using J dsD(e) — 1 we find the implicit expression for 
the potential 



U{B) = Tr lnL4 + B/u(3] - Tr / deD(e) \n(A 
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that has to be evaluated at the saddle point with respect 
to A. 

Liquid state. In the liquid phase we take the replica 
symmetric ansatz A a $ = aS a ^p,B = bS a ,/3, exclude the 
chemical potential p and get 

j3F = [L(a) - ln(a + b/u/3) + bp{\ - p)]/2 - S , (9) 

where L(a) = f dsu(e) ln(a + e 2 ) and So = —p\np— (1 — 
p) ln(l — p) is the high temperature entropy. The corre- 
sponding energy is E — Tp(l — p)b/2. Numerical solution 
of the saddle point equations show that at any density 
the entropy becomes negative at low temperatures. The 
dependence of the temperature at which the entropy be- 
comes zero Ts=o on the density for the coupling constant 
u = 1 is shown in Fig. 1. This line provides the lower 
bound for the glass transition temperature. 

Glass state. In the glass phase the replica symme- 
try is broken, we assume that the solution has one step 
RSB structure and then verify that it is indeed a sta- 
ble solution. Taking B a p = b\5 a $ — 62 Ra,p, A*,/3 — 
Gi<5a,/3 — a2 R a ,p where the matrix R is a block-diagonal 
matrix consisting of m x m blocks with all elements equal 
1, we get the free energy functional 

„„ m — 1 . . , b\ .. L(a\ — ma-A 

PF repl = -^[L{a x) - ln( fll + -)] + - 

- — - XnX \nZ m + pp+-p 2 (b 1 - mb 2 ) (10) 
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where X = a± — ma + (61 — mb^jufi and 
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FIG. 1. Dynamical (solid line) and thermodynamical 
(dashed line) transition temperatures as functions of density. 
The dotted line shows the temperature at which the entropy 
of the liquid state becomes negative. The insert shows the 
trajectory T(m) for thermodynamic (solid line) and marginal 
(dashed line) solutions. The dotted lines represent the tra- 
jectories with constant configuration entropy: The upper one 
(Sconf — 0.07286) corresponds to the physical cooling pro- 
cess in which the system is trapped at the dynamical transi- 
tion temperature. The lower one (S CO nf = 0.04) represents a 
process when a system trapped in a lower energy state. 

with h = p + p(bi — 77162) — b\j2. The free energy ( |10| ) 
has to be taken at the saddle point with respect to 
ai,a,2,bi,b2. For the thermodynamic solution it should 
be also at the saddle point with respect to tti, according 
to (nl) it means that S con f = 0. The numerical solution 
for the dependence of temperature on 777 at fixed den- 
sity (p — 0.5) is shown in Fig 1 insert. By definition 
< 777 < 1, thus the value of temperature for which 
777 = 1 determines the thermodynamic glass transition 
temperature T c . It is plotted as a function of density in 
Fig. 1. 

Stability and marginal solution. In order to analyze 
stability of one step RSB ansatz we expand the Eq.(|^) 
to the second order in fluctuation of the order parameter 
SB and consider different families of fluctuation matri- 
ces SB. This calculation is very similar to the analysis 
of the stability of paramagnetic solution and Parisi solu- 
tion in SK model JlllJl^] so we only sketch it here. We 
find that the most dangerous direction in the fluctuation 
space corresponds to the "replicon" modes fllJlL^] that 
are fluctuations within diagonal blocks of SB satisfying 
the conditions (SB R) a .p = 0, 8B a ^ a — 0. The eigenvalue 
corresponding to these modes is 



A = 



— + 2(52-2.93 + 9i) 



) m /V2n, 



(u/3^ X - l (ai)-ip-92)- 2 
where gk — f P m (z)[l + cxp(—h — z\fb~2)\~ k dz j Z m and 



x(a) — Jexp[— e 2 /2](a + e 2 )~ 2 dx/y/2n. Numerical so- 
lution shows that for any density the thermodynamic 
solution is stable in the wide temperature range below 
T c but eventually becomes unstable at low temperature. 
The solution corresponding to marginally stable states is 
obtained by taking A = 0, the corresponding T(m) de- 
pendence is plotted in Fig. 1 insert. As was explained 
above m = 1 defines the dynamical glass transition tem- 
perature T g , which is plotted as a function of density in 
Fig. 1. 

Physical cooling. During the physical cooling process 
the system remains trapped in a single metastable state. 
As was explained above, the physical properties of such 
system are described by the trajectory m(T) satisfying 
S C onf(iTi,T) = const. We show two such trajectories in 
the insert of Fig. 1. The upper dotted line correspond to 
a typical process in which the system got trapped in the 
highest energy state at the dynamical transition temper- 
ature (for this curve S con f — 0.072864), the lower one to 
a system trapped in a lower energy state corresponding 
to S con f — 0.04. We see that these trajectories cross the 
marginal stability line at lower temperatures, beyond this 
line the one step RSB ansatz becomes unstable. Physi- 
cally, it means that a metastable state in which the sys- 
tem is trapped starts to divide into new ones. 

Thermodynamics of the glass transition. General ther- 
modynamic arguments show that static properties (spe- 
cific heat, pressure) of the glass transition are controlled 
by the total entropy Stot = S + Sconf- We show the 
temperature dependence of Stot and entropy of the liq- 
uid state for density p — 1/2 in Fig 2. The upper dot- 
ted line represents the physical cooling that begins from 
the dynamical glass transition temperature T g . We see 
that the entropy has no jump at the phase transition 
while its derivative (specific heat) does. We also con- 
sider the pressure defined, as usually, by P = —BF/dV. 
Along the physical process S con f = const, therefore 
P = —dF rep /dV. The corresponding plots are shown in 
Fig. 2. As in case with the entropy, the pressure has no 
jump at the phase transition, while its derivative dP/dT 
does. 

In conclusion, we formulated a simple model of the 
structural glass and solved it in the limit of d — > oo using 
the mean field theory. We constructed a full thermo- 
dynamic description of the physical cooling process in 
which the system is trapped in a single metastable state. 
At very low temperature a metastable state in which the 
system is trapped starts to divide into new states, be- 
low this temperature our method based on one step RSB 
ansatz can not be applied. From our analysis we conclude 
that in this system there is no jump in the entropy or 
pressure at the glass transition temperature while there 
are jumps in their temperature derivatives. These qual- 
itative conclusions are in agreement with an established 
phenomenological picture of the glass transition p3| . 
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FIG. 2. Temperature dependence of the entropy and pres- 
sure in the glass phase (solid and dotted lines) and in the 
liquid phase (dashed lines). Solid lines correspond to the 
thermodynamic solution (S CO nf = 0), the upper dotted lines 
represent the total entropy Stot and pressure of the physical 
cooling process (S con f — 0.072864) in which the system is 
trapped at the dynamical transition temperature, the lower 
dotted lines correspond to a process when a system trapped 
in a lower energy state (corresponding to S con f — 0.04). 
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